An integrated procedure to evaluate hydrological models

The growing concern for health‐related problems deriving from pollutants leaching is driving national and international administrations to support the development of tools for evaluating the effects of alternate management scenarios and identifying vulnerable areas. Cropping systems models are powerful tools for evaluating leachates under different environmental, social, and management conditions. As percolating water is the transport vehicle for pollutants transport in soil, a reliable evaluation of water balance models is a fundamental prerequisite for investigating pesticides and nitrate fate. As specific approaches for the evaluation of multi‐layer evolution of state variables are missing, we propose a fuzzy‐based, integrated indicator (I_SWC: 0, best; 1, worst) for a comprehensive evaluation of soil water content (SWC) simulations. We aggregated error metrics with others quantifying the homogeneity of errors across different soil layers, the capability of models to reproduce complex dynamics function of both time and soil depth, and model complexity. We tested I_SWC on a sample dataset where the models CropSyst and CERES‐Wheat were used to simulate SWC for winter wheat systems. I_SWC revealed that, in the explored conditions, the global assessment of the two models' performances allowed identification of CropSyst as the best (average I_SWC = 0·441, with a value of 0·537 obtained by CERES‐Wheat), although each model prevailed for some of the metrics. CropSyst presented the highest accuracy (average agreement module = 0·400), whereas CERES‐Wheat's accuracy was slightly worse, although achieved with a simplified modelling approach (average Akaike Information Criterion = − 230·44), thereby favouring large‐area applicability. The non‐univocal scores achieved by the models for the different metrics support the use of multi‐metric evaluation approaches for quantifying the different aspects of water balance model performances. Copyright © 2010 John Wiley & Sons, Ltd.     

平均经圈环流在大气角动量和感热收支中的作用

本文研究相对座标系中近地表角动量的制造和自由大气中角动量的辐散之间的平衡机制,并与绝对参考系中的机制进行比较。 对1982年元月全球角动量和感热的收支分析证明,被激发的平均经圈环流具有下述特征:其下部水平支施加于大气的惯性矩与近地表角动量的制造相平衡,上部水平支反号的惯性矩与涡动角动量通量的辐合(散)相平衡;而其垂直支的绝热加热效应则平衡着自由大气中热通量的辐散和/或大气的非绝热加热。因而,无论动力上还是热力上,平均经圈环流在大气总环流中都起着重要的作用。     

Potential merits for substorm research from imaging of charge-exchange neutral atoms

The in situ observations of the Earth magnetosphere performed over the past decades of space research have provided a rather good understanding of many partial localized processes of the magnetospheric substorm. The continuing lack of global observations inhibits the construction of a coherent picture of the substorm as a whole, which is actually determined by the coupling of the partial processes. In this context the importance of global observations for the advancement of magnetospheric substorm studies is critical. This paper presents briefly a promising technique of global observations, namely the imaging of charge exchange neutral atoms, or neutral atom imaging (NAI) of the magnetosphere. Model and theoretical estimates of charge-exchange neutral atom fluxes, as well as appropriate spacecraft orbit and instrumentation requirements are presented and discussed for specific regions of interest and vantage points. The potential merits of NAI for substorm research are presented along with possible combinations with other types of observational methods. Substorm issues that would benefit from NAI should include among others the assessment of the ionospheric contribution to the hot magnetospheric plasma, the relative importance of various ionospheric ion source regions, the resolution of spatial and temporal characteristics of substorm ion injections. NAI observations can be precious complements to local observations and lead to the understanding of how local processes, many of which are resolved quite well today, combine to form the global process of the magnetospheric substorm.     

On estimating frequency-magnitude relations from heterogeneous catalogs

The paper studies the effect of magnitude errors on heterogeneous catalogs, by applying the apparent magnitude theory (seeTinti andMulargia, 1985a), which proves to be the most natural and rigorous approach to the problem. Heterogeneities in seismic catalogs are due to a number of various sources and affect both instrumental as well as noninstrumental earthquake compilations.The most frequent basis of heterogeneity is certainly that the recent instrumental records are to be combined with the historic and prehistoric event listings to secure a time coverage, considerably longer than the recurrence time of the major earthquakes. Therefore the case which attracts the greatest attention in the present analysis is that of a catalog consisting of a subset of higher quality data, generallyS ₁, spanning the interval T ₁ (the instrumental catalog), and of a second subset of more uncertain magnitude determination, generallyS ₂, covering a vastly longer interval T ₂ (the historic and/or the geologic catalog). The magnitude threshold of the subcatalogS ₁ is supposedly smaller than that ofS ₂, which, as we will see, is one of the major causes of discrepancy between the apparent magnitude and the true magnitude distributions. We will further suppose that true magnitude occurrences conform to theGutenberg-Richter (GR) law, because the assumption simplified the analysis without reducing the relevancy of our findings.The main results are: 1) the apparent occurrence rate exceeds the true occurrence rate from a certain magnitude onward, saym _GR; 2) the apparent occurrence rate shows two distinct GR regimes separated by an intermediate transition region. The offset between the two regimes is the essential outcome ofS ₁ being heterogeneous with respect toS ₂. The most important consequences of this study are that: 1) it provides a basis to infer the parameters of the true magnitude distribution, by correcting the bias deriving from heterogeneous magnitude errors; 2) it demonstrates that the double GR decay, that several authors have taken as the incontestable proof of the failure of the GR law and of the experimental evidence of the characteristic earthquake theory, is instead perfectly consistent with a GR-type seismicity.     

Evaluation of the SEDD model for predicting sediment yield at the Sicilian experimental SPA2 basin

In this paper a spatially distributed model of the hillslope sediment delivery processes, named the sediment delivery distributed (SEDD) model, is initially reviewed; the model takes into account the sediment delivery processes due to both the hillslope sediment transport and the effects of slope curvature. Then the rainfall and sediment yield events measured at the experimental SPA2 basin, in Sicily, are used both to calibrate the SEDD model and to verify the predictive capability of the distributed sediment delivery approach at event scale. For the SPA2 basin discretized into morphological units and stream tubes, the SEDD model is calibrated at event scale using the measurements carried out at the outlet of the experimental basin in the period December 2000–January 2001. The model calibration is used to determine a relationship useful for estimating the unique coefficient β_e of the model by rainfall erosivity factor R_e at event scale. To test the predictive capability of the β_e = f(R_e) relationship, 20 events measured in the period September 2002–December 2005 are used; the comparison between measured sediment yield values and calculated ones for all monitored events shows that the sediment delivery distributed approach has a good predictive ability at event scale. The analysis also shows that estimating β_e by the relationship β_e = f(R_e) gives a better agreement between measured and calculated sediment yields than obtained with the median value β_e,m of all 27 calculated β_e values. Finally the analysis at annual scale, for the period 2000–2005, allows the estimation of the median value β_a,m representative of the annual behaviour. This analysis shows that the sediment delivery distributed approach also has a good predictive ability at annual scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Effectiveness of simple approaches in mitigating residual deformations in buildings

Developments in performance‐based seismic design and assessment approaches have emphasized the importance of considering residual deformations. Recent investigations have also led to a proposed direct displacement‐based design (DDBD) approach which includes an explicit consideration of the expected residual deformations as an integral part of the design process. Having estimated the expected residual deformations in a structure, engineers are faced with the problem of reducing them to meet the targeted performance levels under pre‐defined seismic hazard levels. Previous studies have identified the post‐yield stiffness as a primary factor influencing the magnitude of residual deformations in single degree of freedom and multiple degree of freedom structures. In this paper, a series of simple approaches to increase the post‐yield stiffness of traditional framed and braced systems for the purpose of reducing residual deformations are investigated. These methods do not utilize recentring post‐tensioned technology. This contribution addresses the feasibility of altering the lateral post‐yield stiffness of structural systems by: (i) using different reinforcement materials with beneficial features in their stress–strain behaviour; (ii) re‐designing the section geometry and properties of primary seismic‐resisting elements; and (iii) introducing a secondary elastic frame to act in parallel with the primary system. The efficiency of each of these techniques is investigated through monotonic and cyclic moment‐curvature and non‐linear time‐history analyses. Of these approaches the design and introduction of an elastic secondary system was found to be most effective and consistent in reducing residual deformations. A simplified design approach for achieving the desired increase of a system's post‐yield stiffness is also presented. Copyright © 2007 John Wiley & Sons, Ltd.     

The June 4th 1999 severe weather episode in San Quirino, Italy: a tornado event?

On the morning of June 4th 1999, a severe weather event took place in San Quirino, a small village of Friuli-Venezia Giulia in the northeast of Italy. This village is located near the piedmont of the Alps, 40 km west from Udine and 60 km north from Venice.Around 0900 UTC (1100 local time), a thunderstorm with an intense hail fall affected the area of San Quirino. A few minutes later (around 0920 UTC, source: a farmer), a funnel cloud from a cumulonimbus touched the ground, producing damages to houses, trees and sheds. The damaged area was quite narrow (about 300 m) and short (less than 10 km). No injuries to people were reported.In spite of the smallness of the area interested by the phenomenon, this storm is studied here starting from the synoptic scale, moving to the mesoscale and finishing with the storm scale, trying to underline its characteristics. These analyses, especially those coming from the Doppler radar images, bring us to the conclusion that the San Quirino episode was produced by a supercell storm.     

Modelling of ground motion in the vicinity of massive structures

A two-dimensional elastic Chebyshev spectral element method (SPEM) is used to model the seismic wavefield within a massive structure and in its vicinity. We consider 2-D models where a linear elastic structure, with quadrangular cross-section, resting on an elastic homogeneous half-space, is impinged upon by the waves generated by a surface impulse at some distance. The scattering of Rayleigh waves and the response of the structure are extensively analysed in a parametric way, varying size, mechanical parameters and shape of the load. Some of the models considered are representative of embankments and earth dams. The simulation shows that some models resonate, storing part of the incoming energy. With realistic parameters, the lowest resonance frequency is due to pure shear deformation and is controlled by the shear velocity and height of the load. Flexural modes are excited only at higher frequencies. The acceleration at the top of the structure may be five/seven times higher than at the base, depending on the mass of the structure. The gradual release of trapped energy produces a ground roll lasting several seconds after the wave front has passed. The ground-roll amplitude depends on the sturcture's mass and can be as large as 30% of the peak acceleration. Outside resonance conditions, the ground motion is almost unaffected by the presence of the artefact; the horizontal motion on top of it is nearly twice the motion at ground level. Similar results should be expected when the incident field is an upcoming shear wave. A qualitative discussion shows that the presence of anelastic attenuation in the embankment does not significantly alter the preceding conclusions, unless it is of very low values (e.g. Q < 15).The modelling results that we discuss indicate that the soil-structure interaction may substantially alter the ‘free-field’ ground motion. From a practical point of view, the main conclusions are: (1) careful analysis is necessary when interpreting seismic records collected in the vicinity of large artefacts; (2) seismic hazard at a site may depend on the presence of man-made structures such as embankments, dams, tall and massive buildings.     

A viscoelastic shear zone model of compressional and extensional plate boundaries

A model is proposed describing the mechanical evolution of a shear zone along compressional and extensional plate boundaries, subject to constant strain rate. The shear zones are assumed as viscoelastic with Maxwell rheology and with elastic and rheological parameters depending on temperature and petrology. Stress and strain are computed as functions of time and depth. For both kinds of boundaries the model reproduces the existence of a shallow seismogenic zone, characterized by a stress concentration. The thickness of the seismogenic layer is evaluated considering the variations of shear stress and frictional strength on faults embedded in the shear zone. Assuming that a fault dislocation takes place, the brittle-ductile transition is assumed to occur at the depth at which the time derivative of total shear stress changes from positive to negative values. The effects of different strain rates and geothermal gradients on the depth of the brittle-ductile transition are studied. The model predictions are consistent with values inferred from seismicity data of different boundary zones.